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ABSTRACT 

We present new observations of the FUV (1100-2200 A) radiation field and the near- to mid-IR 
(3-13.5 //m) spectral energy distribution (SED) of a sample of T Tauri stars selected on the basis of 
bright molecular disks (GM Aur, DM Tau, LkCal5). In each source we find evidence for Ly a induced 
II2 fluorescence and an additional source of FUV continuum emission below 1700 A. Comparison of 
the FUV spectra to a model of II2 excitation suggests that the strong continuum emission is due to 
electron impact excitation of II2. The ultimate source of this excitation is likely X-ray irradiation 
which creates hot photo-electrons mixed in the molecular layer. Analysis of the SED of each object 
finds the presence of inner disk gaps with sizes of a few AU in each of these young (~1 Myr) stellar 
systems. We propose that the presence of strong H2 continuum emission and inner disk clearing are 
related by the increased penetration power of high energy photons in gas rich regions with low grain 
opacity. 

Subject headings: accretion, accretion disks — astrobiology — astrochemistry — circumstellar matter — 
stars: pre-main sequence — ultraviolet :stars 
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1. INTRODUCTION 

In recent years there has been growing evi- 
dence for evolution of solid particles in young 
(< 10 Myr) proto-planeta ry accretion disks 
teeckwith. Henning. fc Nakaeawal l2000t iD'Alessiol 
120031 and references therein). The onset of this evolu- 
tion lies in the coagulation of sub-micron sized particles 
into larger grains; a process which continues until the 
larger grains decouple from the gas and settle to a dusty 
mid-plane. In the mid-plane these solid particles grow 
in size until they become large enough to gravitationally 
focus collisions with smaller bodies, eventually making 
planets (Safronov 1972; WeidenschiUing 1997). 

Since dust grains within the disk absorb stellar UV and 
optical photons and re-emit at wavelengths longer than 
> 2/im, the dust evolutionary process has direct con- 
sequences on disk continuum emission. Thus, at early 
evolutionary stages the presence of dust in the inner disk 
(< 10 AU) is revealed by optically thick emission at near 
and mid-infrared (IR) wavelengths. As grains grow the 
opacity at these wavelengths decreases, revealing stellar 
photospheric emission. The formation of giant planets 
can also affect the dust emission. Gravitational inter- 
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action between the disk and the forming planet results 
in the formation of a gap as the mass of the planet in- 
creases (Bryden et al. 1999; Alibert, Mordasini & Benz 
2004 and references therein), producing a significant de- 
crease in disk flux in the near-IR (Rice et al. 2003, R03). 

What is less recognized, at least in terms of an ob- 
servational signature, is that the molecular evolution is 
closely linked with the dust evolution. Disks in a more 
advanced degree of dust evolution will be more easily 
permeated by the destructive short-wavelength radiation 
fields generated in part by accretion. Since grain evolu- 
tion and planet formation proceeds more rapidly in the 
denser inner disk (WeidenschiUing 1997), these effects 
will be magnified in the very regions that are closest to 
the source(s) of radiation. In this Letter we suggest that 
a strong H2 UV emission feature is an observational con- 
sequence of grain growth/planet formation in the inner 
disks of young 10^ yr) accreting T Tauri disks. 

2. OBSERVATIONS 

The sources chosen for the STIS UV study are 
DM Tau, GM Aur, and LkCa 15. In Table 1 we pro- 
vide some basic characteristics for these systems. All 
sources have relatively similar properties and were se- 
lected on the basis of the presence of gas disks with rich 
molecular complexity (Koerner. Sarecnt. & Beckwith] 
[1293; Dutrey, Guilloteau. & GueUn 1997; Oi et al. 2003}. 
Each are sin gle star systems with n o evidence for binary 
companions Iwhite Gh'^l2nol . HST/STIS spectra 
of DM Tau, LkCa 15, and GM Aur were obtained for 
HST program G09374 on Feb. 5, Feb. 13, and Apr. 
1, 2003 respectively. Exposures were taken using the 
G140L (1150 A - 1730 A) and G230L (1570 A - 3180 
A) gratings with an aperture size of 2". The spectral 
resolution per pixel for G140L is AA = 0.6 A and AA = 
1.58 A for G230L. With a FWHM of the PSF of 1.5 pix 
at the FUV and 2 pix at the NUV, this results in ef- 
fective resolutions of 0.9 A (R~1550) at the FUV and 3 
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A in the NUV (R~ 770). The plate scale is 0.024" per 
pixel in both FUV and NUV, so the aperture size is large 
enough for spectrophotometry. Exposure times were 10m 
(G230L) and ~65m (G140L) for DM Tau and GM Aur 
and 189m (G140L) and 37m (G230L) for LkCa 15. Stan- 
dard CALSTIS pipeline procedures were used to reduce 
the data.^ 

The targets were observed in the mid-IR over the 
course of three separate observing runs (Jan. & Feb. 
2003) that overlapped the HST observations for DM Tau 
and LkCa 15 (5-6 week difference for GM Aur). Ob- 
servations were obtained with the Aerospace Corpora- 
tion's Broadband Array Spectrograph System (BASS; 
Sitko, Lynch, & Russell 2000) with a 3.4" beam on the 
NASA Infrared Telescope Facility. This instrument uses 
a cold beam-splitter to separate the light into two sepa- 
rate wavelength regimes (2.9-6 /im; 6-13.5 /im). Each 
beam is dispersed onto a 58-element Blocked Impu- 
rity Band linear array, thus allowing for simultaneous 
coverage from 2.9-13.5 /im. The spectral resolution is 
wavelength-dependent, ranging from R ^ 30 to 125 over 
each wavelength region. Integration times were 40 min 
(GM Aur), 227 min (DM Tau), and 143 min (LkCal5). 
All observations are calibrated relative to a Tau, with 
typical airmasses ~ 1.07 and calibrator airmasses rang- 
ing from 1.03-1.12. 

3. RESULTS 
3.1. UV Spectra of T Tauri Stars 

Figure 1 shows the FUV de-reddened fluxes of the 
T Tauri stars in our sampl e inclu ding the spectra of 
TW Hya from lHerczeg et all ()2002D . In Table 1 we pro- 
vide the strength of the FUV radiation at 100 AU (esti- 
mated by integrating the FUV radiation field from 1100 
to 1700 A) normalized to that of the interstellar UV ra- 
diation field. This estimate is a lower limit due to the 
unknown Ly a flux.^ 

Figure 2 shows the spectrum of LkCal5 with promi- 
nent emission features identified. In TW Hya (Figure 1) 
there is a strong Ly a emission line that is not as evident 
in the other sources, because these sources are embedded 
within molecular clouds; in contrast TW Hya has no lo- 
cal cl oud and little interstellar absorption ijHerczeg et al.l 
l2004j) . Evidence for strong Ly a radiation that pene- 
trates to the molecular layer can be inferred from the 
presence of numerous emission features coincident with 
Ly a pumped H2 emission lines. In Figure 2 some of 
these coincidences are denoted for LkCa 15, but simi- 
lar features are seen in all sources shown in our sample. 
Beyond the clear H2 emission line features there also ap- 
pears a sharp rise in emission below 1700 A. This feature 
is likely due to a combination of discrete and continuum 
emission emitted by H2 molecules in excited electronic 
states. 

^ There is some overlap in wavelength coverage between the NUV 
and FUV detectors. No correction was needed to be applied to 
match flux levels, and the data is presented with overlap included. 
Some data was not shown at the short wavelength end of the NUV 
spectra to better illustrate particular features. However, the NUV 
data in this spectral region closely mirrors the FUV data. 

^ There is little information of the strength of the radiation field 
below 1100 A. FUSE data of TW Hya indicates that the FUV field 
strength. Go = 3400 (includi ng Ly a), does not ch ange appreciably 
with this radiation included IHerczeg et al.ir2004l) . 
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Fig. 1. — HST/STIS spectra of object s in our sample, inclu ding 
the spectra of TW Hya first discussed in lHerczeg et al] 12003) . but 
with reduced spectral resolution (1.5 A). 



H2 emission below 1700 A can result from at least two 
physical mechanisms (cf. Herczeg et al. 2004), Ly a 
induced fluorescence and electron impact excitation fol- 
lowed by fluorescence (Liu et al. 2002). Ly a exci- 
tation involves only the ungerade electronic states (B 
and C ^H„), producing a spectrum permeated by 
discrete emission lines with some continuum emission. 
However, electron impact excitation can also excite the 
gerade states (E, F etc.) which cascade toward the B, 
C states and emit subsequent UV fluorescence. This 
produces a broader spectrum with greater contribution 
from the H2 dissociation continuum (Abgrall et al. 1997; 
Jonin et al. 2000; Liu et al. 2003). A model spectrum for 
lOOcV electron impact spectrum is displayed in Figure 2. 
A simple comparison of our data to the electron impact 
model in Figure 2 shows similarities that are suggestive 
that electron impact excitation may be important.^ The 
implications of this result are discussed in §4. 

3.2. Spectral Energy Distributions 

The targets, selected on the basis of strong disk molec- 
ular emission, show indications of significant dust evolu- 
tion in their inner disks. For example, the K — L colors of 
LkCa 15 and GM Aur are bluer than 86% of the Taurus 
sample (consisting of 51 CTTS with near and mid-IR ob- 
servations from Kenyon & Hartmann 1995, KH95), and 
the K — L color of DM Tau is bluer than 70% of the sam- 
ple. The BASS observations of the targets strengthen 
these indications. In Figure 3, we compare the mid-IR, 
optical, and IRAS data of the targets to the median SED 
of Classical T Tauri stars in Taurus from D'Alessio et al. 

In the electron impact model there is an additional strong fea- 
ture near 1200 A which we do not discuss due to the close proximity 
of Ly a. In a future publication we will compare the FUV spectra 
with H2 excitation models in greater detail. 
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Fig. 2. — Comparison of observed LkCalS spectrum (solid black 
line), with strong features identified, the electron impact model 
spectrum (dotted blue line), and model of Ly a fluorescence (solid 
red line). The model spectra were produced using an excitation 
model and molecular transition probabilities in the discrete and 
continuum range (Abgrall, RouefT, & Drira 1994, 2000). For ease 
in comparison, the fluxes from both observation and models are 
normalized to the flux at 1425 A. 



(1999). The SEDs of the targets show a clear flux deflcit 
in the near-IR relative to the median in Taurus, analo- 
gous to that of TW Hya (of. Figure 3), for which inner 
disk clearing due to the action of a planet opening a gap 
in the disk has been suggested (§1; C02; R03). 

Attempting to get insight into the structure of the in- 
ner region, we use a simple model, following C02 and 
Uchida et al. (2004). The model consists of: (1) A verti- 
cal "wall" in the inner edge of the outer disk, represent- 
ing the far edge of the gap. This wall, located at radius 
Rt with height is frontally illuminated by the star 
and emits as a black body with temperature T^, where 
T„ = Teff/R]'^/2^/'^ (cf. C02). We are ignoring effects 
of geometry of the gap, inclination, and occultation.^ (2) 
An optically thin inner disk, extending from the dust de- 
struction radius to i?f , with optical depth Tmm at 10 /im. 
This region has a mixture of small (~ 0.1 /im) and large 
(~ 2/xm) grains, and we have added organics and troilite 
to the mostly amorphous silicate grains of Uchida et al. 
(2004). The optically thin dust temperatures are cal- 
culated from radiative equilibrium at each radius. The 
results of this simple modeling procedure are shown in 
Figure 3, where we show the spectrum of each individ- 
ual component; model parameters are given in Table 1. 
It is clear, in the case of LkCa 15, the two components 
alone cannot account for the excess emission above pho- 
tospheric fluxes in the near IR; some excess can be seen 
in the case of GM Aur and DM Tau as well. We find 
that this excess can be explained by black body emission 
at Td = 1400K, the dust destruction temperature (see 
Figure 3). We suggest that the innermost region of the 
inner disk is optically thick and has a rim at the dust 
destruction radius from where this emission arises, as is 
the case in other young disks (Dullemond, Dominik, & 

Emission from the wall in the disk is inclination de- 
pende nt. However, the targets have inclinations ~ 30 - 
60° ISimon. Dutrev. fc Guillotcau 2000; Qi et al. 2003), for 
which the wall emission is with in 50% of the maximum 
JDuUemond. Dominik. &: Nattall200ll) . 



Natta 2001; MuzeroUe et al. 2003, M03). 

Several points can be extracted from our analysis: (1^ 
all disks have cleared their inner regions up to few AU.'^ 
(2) The height of the wall Zw at Rt puts an upper limit 
to the height of the outer disk at that radius because the 
wall could be "puffed up" with the enhanced radiative 
heating. Comparing with predictions from disk struc- 
ture calculations, the low values of z^/Rt in GM Aur 
and DM Tau are consistent with models with very sig- 
nificant grain growth (cf. Figure 5 in D'Alessio et al. 
2001); the much lower value of LkCa 15 cannot be ex- 
plained by models with well mixed gas and dust; it prob- 
ably requires a significant amount of dust settling. (3) 
The column density of dust required to produce the sil- 
icate feature is ^ Tmin/ KilO^m) ^ 0.001 grcm~^, 
with K(10/im) Icm^/gr. Using standard expressions 
for accretion disks (cf. M03), we can estimate the column 
density of gas in the inner disk from the mass accretion 
rate M. With M ~ 3 x lO^^M© yr^^ as representative 
of our objects (Table 1), we obtain I^gas ~ 20grcm~^ 
at 1 AU, for a = 0.01, T = lOOK, and M = IMq, with 
corresponding dust column of ^ 0.2grcm^^ (using the 
standard dust-to-gas mass ratio), much higher than de- 
tected. As discussed in C02, one possibility is that the 
remaining dust is locked up in larger bodies with low 
near-IR opacities. (4) We can estimate the height Zdust 
of the rim at the dust destruction radius from the solid 
angle required to fit the emission, assuming a cylindrical 
geometry. The radius of this cylinder would be the dust 
destruction radius, which can be calculated from the stel- 
lar and accretion luminosities (M03, Table 1). We find 
Zdust < (Table 1), with H ~ O.li?, which is much 

lower than the height of the rim previously found in thick 
inner disks (Dullemond et al. 2001; M03). We note that 
black body emission alone can explain the near-IR excess 
(Figure 3.) This suggests that the innermost optically 
thick region has a small radial extent and moreover, may 
be in the shadow of the rim (Dullemond et al. 2001). 
Taken together, all these results suggest that the solids 
in the inner disks of the targets have experienced signif- 
icant evolution. 

4. DISCUSSION 

In a small sample of stars selected solely on the basis 
of strong and diverse disk molecular emission we have 
found two striking results (1) the presence of strong UV 
continuum emission from H2 molecules and (2) evidence 
for grain growth and possible inner disk clearing by plan- 
etary bodies. Although our sample is limited we suggest 
that these two results may be related. The process of 
grain growth and inner disk clearing will certainly en- 
hance the penetration power of both UV and X-ray ra- 
diation, perhaps even leading to regions that are opti- 
cally thin to UV radiation with penetration only weakly 
limited by molecular opacity. The direct or improved 
exposure of molecular gas to high energy radiation can 
certainly produce the H2 UV emission features seen in 
our study, which only require a small molecular col- 

^ The location of the outer edge of the gap Rt for GM Aur is 
larger than that inferred by Rice et al. (2003). These authors only 
had broad band colors, and included the silicate emission as emis- 
sion from the wall. With our better observations, the optically thin 
emission can be separated from the wall optically thick emission, 
which is much lower in the ~ 10/xm region. 
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Fig. 3. — Spectral energy distributions of GM Aur, DM Tau, and LkCa 15: BASS data (small open circles), optical and IRAS bands from 
KH95 (large solid circles). Also shown are the photospheric fluxes (solid triangles joint by dotted line), and the contribution of the wall 
(short-dashed lines), the optically thin inner region (dotted line), and the possible rim at the dust destruction radius (long dashed line). We 
show for comparison the median of Taurus (solid rectangles joined by dashed lines, error bars denote the quartiles of the distribution) from 
D'Alessio et al. (1999). The last panel shows the SED of TW Hya, with data and model from Calvet et al. (2002). Photospheric fluxes are 
calculated from standard colors and spectral types from KH95, scaled to the de-reddened J magnitude of each object. Models discussed in 
§3.2 are shown in solid lines. We do not include the contribution of the outer disk, which becomes important at longer wavelengths that 
those considered here (C02). 



TABLE 1 

Stellar Properties and Model results 



Stellar Properties 




GM Aur 


LkCa 15 


DM Tau 


SpT 


K3 


K5 


Ml 




4730 


4350 


3720 


Ay 


0.14 


0.62 


0.9 


L 


0.83 


0.74 


0.25 


R* 


1.35 


1.51 


1.20 


log Lace 


-1.149 


-1.055 


-1.749 


log M 


-8.563 


-8.096 


-9.103 


Go 


340 


1500 


240 




Model R 


3SUltS 




Ri/AV 


6.5 


3 


4 


T„/K 


124 


177 


117 


Zw/R 


0.12 


0.03 


0.13 




0.005 


0.007 


0.007 


Rd/R* 


12 


9.8 


12 


Zdustl Rd 


0.01 


0.08 


0.006 



Note. — Go = Field at 100 AU in Habing, 1 Habing = 
standard UV flux from interstellar radiation fleld. Stellar 
data in solar units from Kenyon & Hartmann (1995). 



umn (< 10^^ cm"^: iHerczeg et all lj2004j) '). Additional 
FUV spectra of young TTS combined with forthcoming 
Spitzer/IRS and BASS data on the mid-IR SED will be 
required to fully determine the relation between FUV 
emission and grain growth. 

Electron impact excitation requires the mixture of hot 
electrons within an undissociated H2 layer (Raymond, 
Blair, & Long 1997). For these systems we can as- 
sume that the upper layers of the inner disks are mostly 
cleared of small grains and H2 self-shielding limits photo- 
destruction. Under normal circumstances radiation will 
erode the H2 layer, but these systems are still accreting, 
which provides a source term. In the shielded layer some 
hot photo-electrons will be produced via UV ionization of 
C+j but these do not have sufficient energy 2 eV) for 
H2 excitation. Thus, electron excitation likely requires 
X-ray photons which have greater penetration power and 
produce higher energy electrons. In the dust poor inner 



disk the opacity for abundant 1 keV X -rays will also be 
reduced by nearly a factor of 3 ijGlassgo ld. Naiita. fc Igeal 
IIQQTT allowing X-rays to penetrate into the self-shielded 
layer, ionize H2, and produce hot electrons mixed within 
molecular gas. The radiation produced by electron im- 
pact excitation of H2 , therefore does not originate from 
the star (with an oblique angle of incidence on the disk 
like Ly a), but rather arises from the upper layers of 
the disk itself. This increases the penetration power of 
the UV radiation with resulting effects on disk physics 
(UV clearing of inner disk, accretion, gas temperature 
structure, etc.) and the chemistry of species sensitive to 
radiation below 1700 A. 

The reduction of dust opacity inside the planet form- 
ing regions of the inner disk along with the presence of 
gas (inferred from accretion) suggests that some of the 
heating mechanisms believed important for the outer disk 
will not be effective in the inner disk. In particular, the 
efficiency of photoelectric heating can be expected to be 
reduced and a greater role must be played by X-ray heat- 
ing and UV excited H2 collisional de-excitation. Thus 
the gas physics of an evolving inner disk will change as 
a function of dust evolution. UV emission from H2 must 
trace these changes and provides a direct tracer of the gas 
conditions within the planet-forming regions of T Tauri 
disks that are only weakly probed by the dust. 
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